Eliashberg theory of excitonic insulating transition in graphene 
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A sufficiently strong Coulomb interaction may open an excitonic fermion gap and thus drive 
a semimetal-insulator transition in graphene. In this paper, we study the Eliashberg theory of 
excitonic transition by coupling the fermion gap equation self-consistently to the equation of vacuum 
polarization function. Including the fermion gap into polarization function increases the effective 
strength of Coulomb interaction because it reduces the screening effects due to the collective particle- 
hole excitations. Although this procedure does not change the critical point, it leads to a significant 
enhancement of the dynamical fermion gap in the excitonic insulating phase. The validity of the 
Eliashberg theory is justified by showing that the vertex corrections are suppressed at large N limit. 
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The low-energy elementary excitations of graphene are known to be massless Dirac fermions with linear dispersion 
1, 2j. When Dirac fermions are strictly massless, the system has a semimetal ground state and its Hamiltonian 
possesses a continuous chiral symmetry. However, when the Coulomb interaction is sufficiently strong, the massless 
fermion (particle) may combine with an anti- fermion (hole) to form a stable excitonic pair [3|,|J|. As a consequence, the 
Dirac fermion acquires a finite mass gap, which dynamically breaks the chiral symmetry, and the system undergoes 
a semimetal-insulator transition. In recent years, this excitonic transition has been investigated by several tools, 
including Dyson- Schwinger gap equation Monte Carlo simulation [8rfi~ll]. and renormalization group [Hj]. In 

most of these research works, it was found that such transition can take place when Coulomb interaction strength 
parameter A is larger than certain critical value A c . The opening of fermion gap can result in important changes in 
the low-energy properties [3, Il3l - fil| . 

Currently, there is still no compelling experimental evidence for the existence of excitonic insulating transition in 
graphene. In general, there are three possible reasons why the excitonic transition has not yet been unambiguously 
observed in experiments. First, the Coulomb interaction may be too weak to induce excitonic pair formation. Second, 
the Coulomb interaction is sufficiently strong in certain graphene materials, but the excitonic transition is suppressed 
by various fluctuations and perturbations. Actually, we have recently examined this possibility and showed [5|, [6j that 
even for sufficiently strong Coulomb interaction the excitonic transition can be destroyed by thermal fluctuations, finite 
doping, disorder scattering, and finite lattice effect. Third, the excitonic transition does happen in some graphene, 
but the fermion gap is too small in magnitude to be clearly resolved by any experimental instruments. At present, 
it is not possible to judge what exactly the reason is. In order to examine the third possibility, it is necessary to 
determine the magnitude of fermion gap precisely. 
^ Most previous research efforts focused on an accurate determination of the critical point A c [3I4TTI. fl6j. However, 

an accurate determination of dynamical fermion gap should be equally important. As discussed above, the predicted 
excitonic insulator is experimentally detectable only when the fermion gap is sufficiently large. On the other hand, 
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the graphene with a large gap will have many technological advantages [16j |. As emphasized by Castro Neto [lfj, if 
the fermion gap is too small, the interests of excitonic insulating transition would be purely academic. Therefore, an 
exact determination of the dynamical fermion gap is very important both experimentally and technologically. 

In most previous gap equation analysis of the excitonic transition, the vacuum polarization function appearing in 
the Coulomb interaction function was calculated using the free propagator of massless Dirac fermion, which amounts 
to assuming the random phase approximation (RPA). Thus the feedback effect of the dynamical fermion gap on the 



Coulomb interaction was simply ignored. The importance of such feedback effect can be seen by making a simple 
qualitative analysis. As shown in previous works the screening of Coulomb interaction due to the collective 

particle-hole excitations can have very important influence on the excitonic transition because it reduces the effective 
interaction strength. Technically, the screening of Coulomb interaction is described by the polarization function n(q). 
At the level of RPA, n(q) oc |q|, so the Coulomb interaction is weakened due to dynamical screening. Once a finite 
fermion gap m is included in the polarization function, we have n(q) oc q 2 /m. This term is much smaller than the 
term cx |q| in the low momentum regime, so the effective Coulomb interaction becomes much stronger. From this 
qualitative analysis, we know that the feedback effect of the dynamical fermion gap may play an important role in 
the gap equation analysis of the excitonic transition. 

In order to compute the dynamical fermion gap more accurately it is crucial to study the equations of fermion 
gap and polarization function in a self-consistent manner. This formalism corresponds to the Eliashberg theory 
17], which was originally developed to describe the unusual properties of conventional superconductors with strong 
electron-phonon coupling [17| . Apart from its remarkable success in studying electron-phonon interaction induced 
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superconductors [171 ]. the Eliashberg theory is also useful in many other condensed matter problems. In particular, 
it is widely adopted when stud yin g the non-Fermi liquid behaviors in some correlated electron systems with singular 
fermion-boson interactions [l8l424| . 

In this paper, we study the Eliashberg theory of the excitonic gap generation in graphene. To make a general 
analysis, we include both finite temperature and finite chemical potential. After solving the self-consistent equations 
of fcrmion gap function and polarization function, we found that the critical point does not change, which is reasonable 
because bifurcation theory ensures that the fermion gap can be safely taken to be zero near the critical point. However, 
the size of dynamical fermion gap is significantly enhanced compared to that obtained using the RPA polarization 
function. This implies that the self-consistent Eliashberg theory plays a crucial role in an accurate computation of 
dynamical fermion gap. The Eliashberg theory is justified only when the vertex corrections are unimportant. Within 
the 1/N expansion, we show that the vertex corrections are suppressed in the large N limit. 

The Hamiltonian for interacting Dirac fermions is 

N f - 
i=i Jr 

i N r 2 

+ i^E y rr; ^(r)7o^(r)^^^(r')7o^(r'). (1) 

As usual, we adopt the four-component spinor field ip to describe Dirac fermion and define the conjugate spinor field 
as V> = i/^7o • The 4x4 7-matrices satisfy the Clifford algebra [1, 0]- The physical flavor of Dirac fermion is N = 2. 
The total Hamiltonian possesses a continuous chiral symmetry if) — > e iej5 4>, where 75 anticommutes with 7^, which 
will be dynamically broken once a nonzero fermion gap is generated. 

The free propagator for massless Dirac fermion is Gq 1 (po, p) = joPa — vfI ■ P- After including the interaction effect, 
it is modified to 

G^{Po,p) = ^oPo7o - v F Ai~i ■ p - m(p ,p), (2) 

where Ao,i(po ; p) is the wave function renormalization and m(po, p) is the fermion gap function. These quantities can 
in principle be obtained by solving the following complete Dyson-Schwinger equation 

±-^ 70 G(k)r V(p-k), (3) 

where Tq is the vertex function and V(q) is the Coulomb interaction function. The bare Coulomb interaction is 
Vo(q) = 2 £ u F | q | 1 wn i cn is clearly long-ranged. After taking into account dynamical screening effect from collective 
particle-hole excitations, the effective Coulomb interaction function is modified to 

with n(g) being the vacuum polarization function. The above gap equation has been studied extensively in recent 
years and dynamical gap generation was found when Coulomb interaction is sufficiently strong 0-0] • In most of these 
treatments, only the leading order of 1/N expansion was kept. In particular, both wave function renormalzation and 
vertex corrections were neglected 043) so that Aq 1 = 1 and r = 70 ■ Moreover, the massless fermion propagator 
is widely used when calculating the polarization function shown in Fig.l and the feedback effect of fermion mass is 
simply ignored. Within this approximation, the polarization function has the form 

N a 2 

The aim of this work is to go beyond the popular RPA calculation and include the dynamical fermion gap back into 
the polarization function shown in Fig.l. The importance of this treatment can be readily seen by making a simple 
qualitative analysis. Generally, the effective strength of Coulomb interaction is characterized by two ingredients: 
interaction parameter A = e 2 N /lQvptQ with eo being the dielectric constant and static/dynamical screening due to 
collective particle-hole excitations. While the former is determined by the substrate of graphene sample, the latter 
is reflected in the polarization function. Within RPA, the polarization function (5) behaves as oc q in the static 
limit go = 0. It vanishes linearly as q — > 0, so the Coulomb interaction remains long-ranged. However, the effective 
Coulomb interaction is weakened by the dynamical screening term cx q, which arises from particle-hole excitations. In 



FIG. 1: The bubble Feynman diagram for the polarization function. The solid line is the free fermion propagator and the wavy 
line is the Coulomb interaction function. 



the chiral symmetry broken phase, the Dirac fermion has a finite mass. Assuming a constant mass m, the polarization 
function becomes 

„, , N 9 / m q 2 — 4m 2 q \ 

no? - -q 2 U-2 + , 3 arcsin rr==o ■ 6 

In the limit q <C W, it is simplified to 

n(q) CX (7) 
m 

At the low energy regime, this contribution is much less than the term oc q of the polarization function obtained 
using the massless fermion propagator (RPA). Physically, this reflects the fact that a finite fermion gap reduces the 
dynamical screening effect. Because of this reduction, the effective Coulomb interaction becomes much stronger. 

In general, the fermion mass gap is not a constant, but depends explicitly on momentum. To examine the feedback 
effect of fermion gap on excitonic transition, we will utilize the Eliashberg formalism and couple the fermion gap 
equation to the equation of polarization function self-consistently. 
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FIG. 2: The polarization operator with the vertex correction, but without self-energy correction to intermediate fermions. 

In the Eliashberg formalism, the vertex corrections are usually ignored. This approximation is well-justified in the 
electron-phonon interacting systems as the Migdal theorem ensures that the vertex corrections are suppressed by a 
small factor m/M, where m is the electron mass and M is the nuclei mass [13] ■ In the present problem, we still have a 
small suppressing factor in the vortex corrections based on the 1 /N expansion. This can be explained by considering 
the vertex correction diagram for polarization shown in Fig. 2. In the large N limit, it is found (details are presented 
in the Appendix) to have the form 
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n v ( 9 ) = — j— * — w n( g ). (8) 



Apparently, we know that 



11(g) JV W 

so the polarization function with vertex correction II v (q) is suppressed by a small factor ha.(N)/N comparing with the 
leading polarization function H(q), and therefore can be neglected in the large N limit. One can verify in a similar 
way that the same is also true with the vertex correction diagram for the fermion self-energy. Similar arguments for 
the suppression of vertex corrections are extensively used in the Eliashberg theories of a plenty of physical problems, 



4 



including fermion-gauge systems [Til . Il9j , quantum critical points in itinerant electron systems [20l - |23j , and electron- 
doped cuprate superconductors [24|. Moreover, the wave function renormalization j4qi(po,p) also contain certain 
powers of 1/iV. Therefore, in this paper we will ignore both vertex corrections and wave function renormalizations. 

In order to make a more general analysis, we consider the problem at finite temperature T and finite chemical 
potential fj, and study how the static screening of Coulomb interaction due to finite T and fi is affected by fermion 
mass. We will work in the Matsubara formalism and replace the fermion energy po by an imaginary frequency, 
Po — > iui n = i(2n + 1)ttT. At finite T, it is convenient to adopt an instantaneous approximation [3, |j] and ignore 
the energy-dependence of the interaction function. At finite /z, the frequency iu> n appearing in fermion propagator 
should be replaced by iuj n — /i. The instantaneous approximation allows us to perform the frequency summation over 
/cq analytically, leading to the gap equation 
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(10) 



where q = p k. We need to compute the polarization function Il(q, T) that contains the dynamical fermion mass 
function m(q). Within the instantaneous approximation, the polarization function is defined as 



N /" r/ 2 k 

n fa) = ~jT.J (2^2 Tr k hoG(co n , k + q) 7o ] , 



where the Dirac fermion propagator G(w„,k) contains dynamical fermion mass m(p). Following the procedure 
presented in a previous work [5|, we can sum over imaginary frequency ui n and then obtain 
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where we defined three parameters 



t = [k 2 + x 2 q 2 + 2xk • q + (1 - x)m 2 (k) 

W(k + q) +x(l -x)q 2 ] 1/2 , 
t' = xq 2 + (2x - l)k ■ q+ (1 - x)m 2 (k) 

+im 2 (k + q) — m(k)m(k + q), 
t" = 2k 2 + xq 2 + (2x+l)k-q+(l-x)m 2 (k) 
+a;m 2 (k + q) + m(k)m(k + q). 

Now we obtain the self-consistent equations (10) and (11) for dynamical fermion mass function and polarization 
function. The dynamical fermion mass can be obtained by solving them numerically. 

Here, we would like to make some remarks on the instantaneous approximation. This approximation was originally 
proposed in the study of dynamical fermion gap generation at finite temperature in QED3 [25l |26| . Technically, when 
calculating the fermion gap equation and polarization function, it is possible to sum over the imaginary frequencies 
in the Matsubara formalism only when the energy (frequency) dependence of polarization function is neglected. 
Otherwise, there will be infinitely many coupled gap equations [25l [26| . Moreover, in order to make a qualitative 
analysis of the screening effects, we need to obtain an (semi) analytical expression for the polarization function, as 
Eq.(ll), which can be derived only within the instantaneous approximation. At zero temperature and zero chemical 
potential, it is formally viable to include the energy dependence of the polarization function. Indeed, we have solved 
the fermion gap equation without assuming the instantaneous approximation [5[, but using the RPA polarization 
function (5). The effects of the dynamical part of the RPA polarization function was also investigated in a recent 
paper of Gamayun et al. Q. Due to the non-relativistic nature of the present Coulomb-interacting system, the 
integrations over energy and momenta have to be performed separately, which substantially increases the time needed 
to perform numerical computations. In the Eliashberg formalism, we need to go beyond the RPA level and solve 
the self-consistently coupled equations for the fermion gap and the polarization function. This requires much more 



1A =0.356 



0.1 0.2 0.3 0.4 0.5 0.6 

1A 



FIG. 3: The A dependence of dynamical fermion gap at zero momentum and zero temperature. At physical flavor N = 2, the 
critical strength is given by 1/A C = 0.356. When A < A c , the excitonic insulating transition can not happen. 
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FIG. 4: Dynamical fermion gap at zero temperature. 



computer time and would significantly reduce the precision of numerical output. In order to retain the necessary 
numerical precision, we use the instantaneous approximation even at zero temperature and zero chemical potential. 

Before performing numerical computation, it is helpful to first qualitatively analyze the effect of finite fermion gap 
on the static screening of Coulomb interaction. As an example, we consider the case of finite chemical potential at 
zero temperature. At T = 0, the fermion gap equation is simplified to 



1 [d*] 
m p = — / — 



k m(k) 6>(\A 2 +m 2 (k) - /i) 



which couples to the polarization function 
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Now we assume a constant fermion mass gap m, then the integration over momentum k in these equations can be 
carried out exactly. When fi > ^Jm 2 + q 2 /4, the polarization function behaves as n(q) = Nfi/n. This expression 
implies that the Coulomb interaction is now statically screened by finite chemical potential and thus becomes short- 
ranged. Such static screening effect will rapidly destroy the excitonic pairing instability However, when the 
fermion gap is relatively large, m > fj,, we have n(q) oc q 2 /m in the low momentum regime, so the Coulomb 
interaction remains long-ranged and is only poorly screened. From this qualitative analysis, we know that once the 
feedback effect of fermion gap is included, the static screening may be suppressed and the effective Coulomb interaction 
may still be strong even at finite chemical potential. Besides chemical potential, thermal fluctuation can also induce 
static screening [5| . The effect of fermion gap on static screening at finite temperature can be qualitatively analyzed 
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similarly. To gain more quantitative understanding on feedback effect of fermion gap, we have to perform numerical 
computation. 
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FIG. 5: Dynamical fermion gap at finite temperature and zero chemical potential, with 1/A = 0.2. 
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FIG. 6: Dynamical fermion gap at finite chemical potential with T/A = 10 7 and 1/A = 0.2. 

We have numerically solved the coupled equations of dynamical fermion gap and polarization function by means 
of straightforward iterative method. When T = 0, the fermion gap of zero momentum is shown in Figj3l where the 
blue (red) curve represents the results obtained using the perturbative (self-consistent) polarization function. From 
the numerical results, we know that the critical Coulomb interaction parameter A c takes the same value in these two 
cases (1/A C = 0.356 for physical flavor N — 2). This fact is easy to understand since the bifurcation theory requires 
that the fermion gap should vanish near the critical point between gapless and gapped phases. However, away from 
the critical point in the insulating phase, the fermion gap obtained using different polarization functions are no longer 
the same. The momentum dependence of fermion gap at T = is shown in FigU for different values of A. The 
solid (dashed) lines are the results obtained using perturbative (self-consistent) polarization function. The blue, red, 
and green lines correspond to 1/A = 0.1,0.2,0.3 respectively. From both Figl3]and FigJH it is easy to see that the 
dynamical fermion gap is significantly enhanced in the Eliashberg formalism. Specifically, for graphene suspended in 
the vacuum with 1/A = 0.295 [l6j], the fermion gap calculated by the Eliashberg theory is greater by two orders of 
magnitude than that calculated using perturbative polarization. 

It is not hard to include the effects of temperature and chemical potential. The dynamical fermion gap at finite 
T and zero /i is shown in FigfSJ and the fermion gap at finite /i with T/A = 10 -7 is shown in FigJFJJ In both 
of these diagrams, the solid (dashed) lines are the results obtained using pertubative (self-consistent) polarization 
function. The interaction strength is fixed at 1/A = 0.2. The blue, red, green, and black lines are the results at 
T/A = 10~ 6 , 5 x 10~ 6 , 10~ 5 , 2 x 10~ 5 respectively in FigE]and arc the results at /x/A = 10~ 6 , 5 x 10~ 6 , 10~ 5 , 2 x 10~ 5 
respectively in FigjB] Here, A is the ultraviolet cutoff for momenta. Apparently, there are substantial enhancement 
effects of dynamical fermion gap in both cases of finite temperature and finite chemical potential. In the Eliashberg 
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theory, the static screening effects caused by finite temperature and finite chemical potential are suppressed, so the 
fermion gap is nearly independent of T (see Fig(5]) and fj, (see FigJS]). 

The fermion gap enhancement can be understood as follows. It is well-known that the screening of Coulomb 
interaction is determined by the density of states of fermions. Once the feedback of dynamical fermion gap is taken 
into account, the density of states of fermions is significantly reduced and the static or dynamical screening effect 
becomes less important. As a consequence, the effective Coulomb interaction becomes stronger, which in turn leads 
to larger fermion gap. 

In summary, we have studied the Eliashberg theory of excitonic phase transition in graphene. After solving the cou- 
pled equations for dynamical fermion gap and polarization function in a self-consistent manner, we found a significant 
enhancement of dynamical fermion gap in the excitonic insulating phase. The enhancement found within Eliashberg 
formalism is owing to the suppression of static or dynamical screening of Coulomb interaction by dynamical fermion 
gap. Therefore, the self-consistent treatment of polarization function should be used in the accurate computation of 
the dynamically generated fermion gap. The validity of the Eliashberg theory is justified by showing that the vertex 
corrections are suppressed by a small factor in the large N limit. 

This work is supported by the National Natural Science Foundation of China under grant No. 11074234. G.Z.L. 
is also supported by the Project Sponsored by the Overseas Academic Training Funds of University of Science and 
Technology of China. 



Appendix 

The polarization function Tl v (q) shown in Fig. 2 is defined as 

f d 3 k f d 3 v 

n v (?) = -N J ^3 J j^Tr[ l0 G a {k + q) la G (k+ P + q)-f a 

xG Q (k+ P ) 7o D( P )Go(k)}. (12) 

To calculate this function, we will follow the method of Franz et al. [27j . We are mainly interested in the leading 
behavior of II v (g) in the q — > limit. In this limit, the above integral has singularities as k — > and k — > —p. Thus, 
we may evaluate the whole integral by expanding the regular parts of the integrand near these two singular points. 
Keeping only the leading terms, we have 

d 3 k f d 3 p 



/d k f dp 



xG (phoD(p)G (k)} 

/d 3 k f d 3 p 

xG (k+p)~f D(p)G (-p)] (13) 
Perform a variable shift, k — > k — p, for the second term, then 



n v (<z) = 2ATr 

where 



d 3 k 

X[<1) I J^ G o( k hoG (k + q) 



d 3 p 
(2n) 



The most leading term is found to be 



,2 , 



In the large N limit, it is possible to use the approximation 

cos 2 6 - sin 2 6 32e Q v F , / e 2 
smf Ne 2 V Se v F 



d0 t Ne2 , HAT— ), (17) 
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*(?) = - / T^TsToGob + qhoGoiphoDip). (15) 
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X(q) = ^-ln - / d9 . (16) 
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so that 



4 70 ^ (N^) In (A) 



(18) 
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After substituting Eq.([T8]) to Eq.([T4]), we finally get 



n v (g) 
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\ 8e v F J \Q J 



U(q). 



(19) 
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